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spread by a plurality of spreading codes, which method 
comprises the steps of: 

(a) converting the CDMA signal to a baseband signal; 

(b) recognizing symbol boundaries in every multipath for 
every user in the baseband signal to produce a delay 
estimate; 

(c) estimating an attenuation and phase rotation impressed 
on the signal by each multipath to produce channel 
estimates; 

(d) generating new symbol estimates for every user at a 
sampling rate that is greater than a symbol rate of the 
signal; 

(e) generating an estimate of the received signal using 
said symbol estimates and channel estimates; 

(f) comparing the received signal with its estimate and 
feeding estimation error back to a means for estimating 
the transmitted symbols; 

(g) sampling the symbol estimates at the estimated sym- 
bol boundaries of each user to obtain final symbol 
estimates; and 

(h) repeating steps (b) through (g) with every new sample. 
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MULTI-USER CODE DIVISION MULTIPLE 
ACCESS RECEIVER 

BACKGROUND OF THE INVENTION 
THIS INVENTION relates to a multi-user code division 
multiple access receiver based on recursive estimation. 
More particularly the division relating to a method of 
receiving a code-division multiple-access signal. 

Code division multiple access (or CDMA) cellular com- 
munications systems suffer from co-channel interference (or 
multi-user or multiple-access interference [MAI]) on the 
reverse or up link (i.e. from a mobile handset to a base 
station). This is due to the impossibility of maintaining 
orthogonality between code channels used by independent 
mobile handsets, which transmit asynchronously. This form 
of interference limits the uplink capacity severely, and very 
significant capacity gains can be achieved if multi-user 
interference can be reduced, or if joint detection of all users 
is employed, as opposed to conventional single-user detec- 
tion techniques currendy in use which merely treat multiple - 
access interference (MAI) as uncorrelated noise. 
. Previous work on MAI cancellation dates back to S. Verd 
u, "Minimum probability of error for asynchronous Gauss- 
ian multiple-access channels," IEEE Trans. Inf. Theory, vol. 
32, no. 1, pp. 85-96, January 1986, which describes the use 
of a Viterbi decoder in concert with a maximum-likelihood 
(ML) decoding metric for the demodulation of the bits 
transmitted by each user. The maximum-likelihood detector 
has a complexity which increases exponentially — in a binary 
phase shift keyed (BPSK) system, this is of the order of 2* 
multiply-and-add's, K being the number of users in the 
system. Its complexity prevents the practical implementa- 
tion of the ML detector, and thus other sub-optimal schemes 
have been proposed. 

Among these, direct descendants of Verdu 's ML detector 
are the sequential decoder discussed in Z. Xie, C. K. 
Rushforth, and R. T Short, "Multi-user signal detection 
using sequential decoding," IEEE Trans. Communications, 
vol. 38, no. 5, pp. 578-583, May 1990, and a limited 
tree-search detector discussed in L. Wei, L. K. Rasmussen, 
and R. Wyrwas, "Near-optimum tree-search detection 
schemes for bit-synchronous multi-user CDMA systems 
over Gaussian and two-path Rayleigh fading channels/' 
IEEE Trans. Communications, Vol. 45, No. 6, pp. 691-700, 
June 1997, which operate on the principle of searching the 
decoding tree structure non-exhaustively, as opposed to the 
ML algorithm which examines every branch of the tree. The 
complexity of such schemes is still fairly large, and their 
structure is more suited to a programmable digital signal 
processor (DSP) implementation than a large-scale inte- 
grated (LSI) circuit one. Present DSP's are still not powerful 
enough to implement the tree-search family of multi-user 
detectors. 

The conventional detector passes the received signal 
through a filter bank, each branch of which is matched to the 
spreading code of one user. The K outputs of the filter bank 
are then passed through K decision devices (or slicers) 
which decide on the symbols that were transmitted. This 
simple technique suffers from MAI as described earlier, but 
an additional matrix filter based on the cross-correlations 
between the spreading codes used by individual users can be 
designed to remove MAI completely. This detector structure 
is known as the decorrelating detector and was proposed in 
R. Lupas and S. Verdu, "Linear multi-user detectors for 
synchronous code -division multiple-access channels/' IEEE 
Trans. Inf. Theory, vol. 35, no. 1, pp. 123-136, January 
1989. 



10,099 Bl 

2 

The problems with the decorrelating detector are that: 

1. when long spreading codes with periods greater than 
the symbol interval are used, the matrix filter has to be 
re-calculated in every symbol. This is a non-trivial 

S operation requiring the inversion of a KxK matrix; 

2. with long codes, there is a non-zero probability that in 
a given symbol interval, the codes selected by the K 
users will not form a linearly independent set. This 
would be fatal to the decorrelating detector; and 

10 3. The decorrelating detector removes MAI without 
regard for noise and interference from other sources, and can 
be shown always to increase the noise component not due to 
MAI. This noise enhancement problem is a characteristic of 
zero-forcing filters such as the decorrelating detector. 

15 The first two problems are almost impossible to overcome 
due to the very concept of the decorrelating detector, which 
is to invert the channel matrix in order to recover the 
transmitted signal vector. 
The noise enhancement problem can be limited using the 

20 same technique as in channel equalisation, namely the use of 
a minimum mean squared error (MMSE) solution. This was 
proposed in Z. Xie, R. T. Short, and C. K. Rushforth, "A 
family of sub-optimum detectors for coherent multi-user 
communications/' IEEE J. Sel. Areas Comms., vol. 8, no. 4, 

25 pp. 683-690, May 1990, and further explored in U. Madhow 
and M. L. Honig, "MMSE interference suppression for 
direct-sequence spread-spectrum CDMA," IEEE Trans. 
Communications, vol. 42, no. 12, pp. 3178-3188, December 

1994, The MMSE detector also in theory requires a matrix 
30 inversion, and in addition to knowledge of spreading codes, 

also requires knowledge of the thermal noise power at the 
receiver. This may seem impractical, but the amount of side 
information required can be reduced to merely a training 
sequence when an adaptive implementation is applied. This 

35 approach is expounded in, for example, S. L. Miller, "An 
adaptive direct-sequence code-division multiple-access 
receiver for multi-user interference rejection," IEEE Trans. 
Communications, vol 43, no. 2/3/4, pp. 1746-1755, 
February/M arch/April 1995, S. Verdu, "Adaptive multi-user 

40 detection," in Code Division Multiple Access 
Communications, S. G. Glisic and P. A. Leppanen, Eds., pp. 
97-116, Kluwer Academic Publishers, The Netherlands, 

1995, and a blind approach not requiring the use of training 
symbols has also been proposed in M. Honig, U. Madhow, 

45 and S. Verdu, "Blind adaptive multi-user detection," IEEE 
Trans. Inf. Theory, vol. 41, no. 4, pp. 944-960, July 1995. 

The adaptive MMSE detector however suffers from one 
major drawback in that it can only be used for short-code 
systems, i.e., those which use spreading codes with periods 

50 exactly equal to one symbol interval. Such systems do not 
currendy exist and are not under consideration for future 
cellular systems because of the need to have codes long 
enough to uniquely identify every single mobile handset in 
use at one time around the world, 

55 Another family of multi-user detectors works on the 
principle of interference cancellation through progressive 
subtraction of estimated interference components from the 
received signal. The successive interference canceller (SIC) 
as discussed in P. Patel and J. Holtzmann, "Analysis of 

60 simple successive interference cancellation scheme in a 
DS/CDMA system," IEEE J. Sel. Areas Comms., vol. 12, 
no. 5, pp. 796-807, June 1994, and the parallel interference 
canceller (PIC) discussed in M. K. Varanasi and B. Aazhang, 
"Multi-stage detection in asynchronous code -division 

65 multiple-access communications," IEEE Trans. 
Communications, vol. 38, no, 4, pp. 509-519, April 1990, 
are well known and many improvements have been made to 
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the basic structures. Notable ones include hybrid serial- (a) converting the CDMA signal to a baseband signal; 

parallel structures as discussed in M. Sawahashi, H. Andob ( D ) recognising symbol boundaries in every multipath for 

and K. Higuchi, "DS-CDMA pilot and data symbol-assisted every user i n th e baseband signal to produce a delay 

coherent multistage interference canceller using repeatedly estimate* 

updated channel estimation," in the Proc. IEEE Int'l Conf. c , x .. . ' „ . . * *• • a 

Comm. Systems CCCS)/Int'l Wkshop Intelligent Sig. Proc. < c > estimating an attenuation and phase rotation impressed 

& Comm. Systems (ISPACS), Singapore, November 1996, on . the ^ b ? each multl P ath to P roduce channel 

pp. 585-589 and the use of tentative decisions at the output estimates; 

of each stage instead of hard decisions as discussed in U.S. (d) generating new symbol estimates for every user at a 

Pat. No. 5,579,304 and D. Divsalar and M. Simon, sampling rate that is greater than a symbol rate of the 

"Improved CDMA performance using parallel interference signal; 

cancellation," Tech. Rep. 95-21, Jet Propulsion Lab., Cali- ( e ) generating an estimate of the received signal using 

fornia Inst, of Tech., October 1995. The linear clipper sa id that is greater than a symbol rate of the signal 

soft-decision function in particular produces quite spectacu- symbol estimates and channel estimates; 

lar improvements in performance over the hard-decision / A mm „.„*„„ r ™™,i o1 ;#e . * 

SIC, I shown in X Zhang and D. Brady, "Soft-decision « (9 comparing the received sipal with its estimate and 

multistage detection for asynchronous AWGN channels", in * es ^ 0D ™ r back to a means for estimaUn S 

the Proc. 31st Annual Allerton Conf. Comms,, Cont. & , , me «™ &miiicd s V mbo ^ 

Comp., Monticello, 111., 1993 and for another algorithm in L. fe) sampling the symbol estimates at the estimated sym- 

B. Nelson and H. V. Poor, "Iterative multi-user receivers for bo1 boundaries of each user to obtain final symbol 

CDMA channels: an EM-based approach," IEEE Trans. 20 estimates; and 

Communications, vol.44, No. 12, pp.1700-1710, December CO repeating steps (b) through (g) with every new sample. 

1996. A code division multiple access signal receiver compris- 

Various forms of the subtractive interference canceller are m g: means to receive a CDMA signal combining symbol 

disclosed in U.S. Pat. No. 5,579,304; 5,218,619; 5,467,368; values spread by a plurality of spreading codes; means to 

and 5,363,403. 25 convert the CDMA signal to a baseband signal; means to 

Due to the large number of different types of multi-user recognise symbol boundaries in every multipath for every 

CDMA detectors which have been proposed, it is quite a task user in the baseband signal to produce a delay estimate; 

to rank them according to performance beyond noting that means to estimate an attenuation and phase rotation 

the ML detector performs best while the conventional detec- impressed on the signal by each multipath to produce 

tor performs worse. Therefore, we do not claim that this 30 channel estimates; means to generate new symbol estimates 

invention performs better than all other multi-user CDMA for every user at a sampling rate greater than a symbol rate 

detectors, only that against a similar detector described in D. of the signal; means to generate an estimate of the received 

S. Chen and S. Roy, "An adaptive multi-user receiver for si S nal ^mg said symbol estimates and channel estimates; 

CDMA systems" IEEE J. Sel. Areas Comms., vol, 12, no 5, means to compare the received signal with its estimate and 

pp. 808-816, June 1994, it does perform a lot better and 35 feed a estimation error back to a means for estimating the 

against a four-stage linear-clipper SIC, performance is com- transmitted symbols; and means to sample the symbol 

parable as will be demonstrated later. estimates at the estimated symbol boundaries of each user to 

In D. S. Chen and S. Roy, "An adaptive multi-user obtain a final symbol estimate, 

receiver for CDMA systems" IEEE J. Sel. Areas Comms., Mobile communication channels change over time, and it 

vol. 12, no 5, pp. 808-816, June 1994, a detector which 40 ^ necessary to track these changes in order to maintain 

operates on the basis of recursive least squares (RLS) weight acceptable receiver performance. A multipath searcher is a 

updating and the concept of regenerating the received signal device which constantly probes the channel to ensure that 

based on known spreading codes and channel coefficients the strongest paths are being used in the receiver, 

was described. The channel model used is similar to the one Traditionally, it consists of one or more sliding correlators, 

described herein, but it assumes synchronous transmission, 45 which correlate the received signal with time-shifted ver- 

i.e., all users' symbol intervals are synchronized, whereas sion s of the spreading code of the desired user. The outputs 

this invention explicitly includes asynchronous systems, of the correlators for aU time delays within a specified range 

Also, the weight up-dating algorithm described herein dif- are recorded and, if the L strongest paths are to be used in 

fers from the RLS algorithm used in D. S. Chen and S. Roy, me receiver, the delays corresponding to the L largest 

"An adaptive multi-user receiver for CDMA systems" IEEE 50 correlator outputs are selected. 

J. Sel. Areas Comms., vol. 12, no 5, pp. 808-816, June 1994. In a multi-user CDMA system, the received signal is 

In T J. Lim and L. K. Rasmussen, "Adaptive symbol and corrupted by MAI, and thus using it directly in a sliding 

parameter estimation in asynchronous multi-user CDMA correlator leads to poorer performance than in a single-user 

detectors", IEEE Trans. Communications, vol. 45, no. 2, pp. system. 

213-220, February 1997, a joint parameter estimation and 55 Accordingly, another aspect of the present invention pro- 
bit detection algorithm using the extended Kalman filter y i dcs a method of multipath searching at the output of a 
(EKF) was proposed. The novelty in that algorithm lay in the receiver, comprising the steps of: respreading the estimated 
joint adaptation of time delay estimates and bit estimates symbols transmitted by each user; summing the respread 
with the EKF, which differs from the Kalman filter used in signals; subtracting the summed signal from the received 
this invention because the latter can be used only when time 60 signal to derive a residual error signal; and feeding the 
delays are assumed to be known. residual error signal to a bank of multipath searchers. 

SUMMARY OF THE INVENTION BRIEF DESCRIPTION OF THE DRAWINGS 

A method for receiving a code-division division multiple- In order that the present invention may be more readily 

access (CDMA) signal combining symbol values spread by 65 understood, embodiments thereof will now be described, by 

a plurality of spreading codes, which method comprises the way of example, with reference to the following drawings, 

steps of: in which: 
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FIG. 1 is a schematic block diagram of a receiver structure cell, h k (i) is the symbol decision for the kth user in the ith 

embodying the present invention; symbol interval or epoch, z is a vector of multipath delay 

FIG. 2 A is a schematic block diagram of another receiver estimates, c is a vector of complex multipath coefficients, T 

structure embodying the present invention; is the symbol interval, z -1 represents a one-sample delay, 

FIG. 2B is a maximal ratio combiner for use with the 5 md s *( n ) te the nth sam P le of me kth user ' s spreading 

receiver structure of FIG. 2A; waveform. 

FIG. 3 is a schematic block diagram of a joint channel Bec * USe indivMual ™" a » assumed to transmit in an 

. a a. a a a. c i * < • . - asynchronous manner, symbol decisions are made at differ- 

estimator and data demodulator for use with embodiments of ' . c a c , \ > i a - i_ 

# . _ . • ent times. For user k^n) is sampled at time n-iT+r*, where 

the present invention „„ ~ . u i j *• j * • /* * \ j ,t_ 

^ 10 T is the symbol duration and T fc =min(T w , . . . , t w ), and the 

FIG. 4 shows bit error rate curves for a receiver structure result fed to a slicer or memory i ess decision device for 

embodying the present invention and a conventional detec- symbol decoding . It should be noted lhat all time variables 

tor in a four-user system with a processing gain of 16; are normalised to the sampling interval. This process is 

FIG. 5 shows bit error rate curves for a receiver structure depicted at the bottom of FIG. 1. 
embodying the present invention and a conventional detec- 15 Another receiver structure embodying the present inven- 
tor in a four-user system with a processing gain of 32; tion is shown in FIG. 2A. For this structure, it is not 

FIG. 6 shows bit error rate curves for a receiver structure necessary to assume that the path delay are all much smaller 

embodying the present invention and for a conventional than T. In this case, in place of the symbol estimate vector 

detector with a processing gain of 16; b(n)=[D 1 (n), . . . , bj^n)] 7 ^ a state estimate vector is updated 

FIG. 7 shows bit error rate curves for an RLS detector, a 20 m '' State Update" block. The state estimate 

conventional detector and a receiver structure embodying vector is: 

the present invention; r w * m* /„\t 

FIG. 8 is a schematic representation of a system embody- 
ing the present invention for multipath searching; and ^ Each element of the state vector is an estimate of one 

FIG. 9 is a schematic representation of a system embody- symbol-modulated path coefficient, e.g. 
ing the present invention for decoding access channels. 



DETAILED DESCRIPTION 



and so forth. In this scheme, channel estimates c Jk are not 
FIG. 1 shows a block diagram of a code division multiple 30 required in the state estimate update algorithm, but are used 
access receiver structure embodying the present invention. f° r coherent demodulation in a maximal ratio combiner, as 
The signal r(n) represents the nth sample of the continuous- shown in FIG. 2B. 

time baseband received signal, sampled at a rate of l/T s . In both the receiver structures of FIGS. 1 and 2A, the key 
From the following expression, an estimate may be obtained point is that a vector is recursively changed or updated in a 
of r(n) and is expressed as f(n): 35 systematic fashion, with the objective of getting progres- 

sively better estimates of the symbols that were transmitted 
k (i) by each user. Such an objective may be met by minimising 



r(n)-- 



Eh(n)Y Zusk(n-tkt) either of the two following functions with respect to the 



vector 6(n) or x(n) (depending on which receiver structure is 

40 used): 

where K denotes the number of active mobile users in the A(»)-^iK")-^(»)li 2 (2) 

cell and L is the number of multipaths per user tracked at the j 2 (n)*E\r(n)-?(n)\ 2 (3) 

receiver. The path delay estimates t w and complex channel 2 

coefficient estimates c kI are assumed to have been obtained 45 In (2), St(n) is understood to include both x(n) in FIG. 2A as 

using an unspecified channel estimator and multipath well as b(n) in FIG. 1. 

searcher/tracker. s^(n) is the spreading code used by the kth The differences between J^n) and J 2 (n) are subtle but 

user in the system and b^n) is an estimate of the symbol significant. Minimising J a (n) requires that information 

transmitted by the kth user in the system which contributes regarding the time variation of x(n) be known, or assumed; 

to the received signal at sample instant n. 50 minimising J 2 (n) requires no such information. However, 

The a priori estimation error e(n)=r(n)-f(n) is fed^ to a this is not to say that the former approach is inferior to the 

recursive update algorithm for b^n), which generates 6^+ latter— in fact the opposite is true because minimising Jj(n) 

1), the next estimate of the transmitted symbol. The updating admits a further degree of freedom in the modelling of the 

algorithm is designed to bring fi^n) ever closer to b^i), the signal, which can be used to advantage as discussed below, 

true transmitted symbol which contributes to the received 55 A large number of adaptive or recursive algorithms may 

signal at the time n. (It is assumed that the largest path delay be derived based on the minimisation of Jj(n) or J 2 (n) in 

is much smaller than the symbol interval T). The index i equations (2) and (3). Examples include stochastic steepest 

denotes the ith symbol interval. In all such algorithms, descent and recursive Gauss-Newton algorithms, already 

parameters which govern its behaviour (such as step size, well known in the literature. 

initial values, etc.) need to be specified, and these are 60 In this example, a Kalman filter-based algorithm is devel- 

identified together under "Parameters" in FIG. 1. Estimates oped for use with the receiver structure of FIG. 1, assuming 

of the multipath delays and complex coefficients may also be that binary phase shift keying (BPSK) modulation is 

needed in the Symbol Estimate Update block, and this is employed. 

therefore indicated in FIG. 1 as well. The Kalman filter is a well-known tool in statistical 

Referring to FIG. 1, all time variables are in units of 65 detection and estimation theory, and has countless variants, 

sampling interval, r(n) is the sampled received signal which have been introduced either to improve numerical 

sequence, S represents the spreading codes of all users in the stability or computational efficiency. The most famous 
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examples are the so-called "square root" Kalman filtering 
algorithms, which update a square -root factor of the cova- 
riance matrix P(n) instead of P(n) itself, in order to preserve 
the symmetry of this matrix. 

Assuming that each of K users transmit over a multipath 
channel with L paths each and that the maximum path delay 
for any user is much smaller than the symbol interval T, the 
received signal is given by: 



(4) 



L 



= S(n)B{n) + v{n) 



where S(n) is a row vector of channel-modified spreading 
codes, fi^n) is defined as 

k (n) = £ b k [()u{n-rr-T,\ 



b^i) is the ith symbol transmitted by user k, u(n) is a unit 
pulse of duration T, x k is the smallest path delay for user k, 
S^(n) is the kth spreading code, c w is the 1th channel 
coefficient for the kth user, the x hl is the 1th path delay in the 
kth channel. v(n) is white Gaussian noise of zero mean and 
variance a 2 . 

It is apparent that 6(n)=b(n-l) except at symbol 
boundaries, i.e., when n=iT+T^jc=l, . . . , K. This time 
variation may be expressed concisely in the equation 

(5) 

where w(n) is equal to the zero vector except when n^iT+Tjt, 
in which case it takes the values 0 and [0, . . . , 0,-25^), 
0, . . . , 0] r with equal probability. This point is easily verified 
by considering that b^i+l^b^i) or-b^i) with equal prob- 
ability. 

The vector w(n) can therefore be seen as a non-stationary 
random vector with a covariance matrix at time n of: 



G(n) = E[w{n)yJ{n)] 



if n * iT r k 



diag(0 0, 7b k (n), 0 0 if n = iT +■ r k 



If it is assumed that w(n) is Gaussian and has zero mean, 
the Kalman filter can be used to iteratively determine the 
MMSE state estimates b(n), using the state-space models (5) 
and (4). The kth element of o(n) sampled at the bit bound- 
aries would yield bit estimates for user k. 

The recursive algorithm based oo this description is 
summarised below. This algorithm requires that the relative 
time delays x w and channel coefficients c w are known, which 
is not an impractical requirement since code synchronisation 
and channel estimation are needed in all conventional 
CDMA receivers anyway. 

Initialisation 

i>(0)-I 

G(o)-o 

fr(0)-0 
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Input parameters 

o^Rough estimate of noise power 
At each update, do: 

?(«)>^n-l)+fi(n-l) 
P{n):-{I-K(n)S(n)]r(n) 

This last section represents the steps needed in the Symbol 
Estimate Update block of FIG. 1. 

Although the present example implementation assumes 
BPSK modulation, complex multilevel modulation schemes 
such as M-ary PSK, quadrature amplitude modulation 
(QAM), differential PSK, etc. can also be employed with 
suitable adjustments to the algorithm. 

The receiver structure shown in FIG. 2 A may also be used 
for multipath channel estimation i.e., the estimation of the 
complex multipath coefficients c w , k=l, . . . , K, 1=1, . . . , 
L, which together with the path delays x kl constitute a 
mathematical model for the physical baseband channel. Data 
demodulation using the well-known Rake combiner as dis- 
cussed in J, G. Proakis, Digital Communications, McGraw- 
Hill, 3rd edition, 1995 and most other more sophisticated 
receivers depends crucially on the availability of accurate 
channel estimates, and thus it is usually an important issue 
in a CDMA receiver. 

To illustrate the concept, FIG. 3 depicts the proposed 
scheme for the kth user assuming transmission over a 
single-path channel, which necessitates the estimation of 
just one channel coefficient ck (t), in general a Rayleigh 
fading complex random variable and so slowly time-varying 
(hence the time index t). x*(n) is the kth element of the state 
estimate vector x(n) generated using the Kalman filter or 
some other recursive algorithm. Sampling this high-rate 
(chip rate or higher) signal at the symbol boundaries i.e., 
when n=iT+x fr where x k is the most recently available time 
delay estimate for this user, yields a signal which should in 
theory be close to the product of the channel coefficient and 
the symbol transmitted, c^ijb^i). In writing this expression, 
it is assumed that c k is constant over one symbol interval, a 
valid assumption provided the maximum Doppler frequency 
is very much smaller than the symbol rate, a condition easily 
satisfied in practice. 

The single-path channel model is easily extended to 
multipath models with known path delays. 

Assuming the use of phase shift keying (PSK) with unit 
energy symbols (in other words the signal constellation 
consists of points which lie on the unit circle in the complex 
plane), multiplying $ k (n) with b^(i)*, where (.)* represents 
complex conjugation, at the symbol boundary a=iY+x k 
should produce an estimate of c^i). Of course, b*(i) is 
unknown except when pilot symbols are being transmitted. 
Outside the pilot-symbol intervals it is proposed to feed the 
estimated b^i), denoted b^i), back in the manner shown in 
FIG. 3, to obtain the channel estimate c fc (i). To further 
improve the performance of the channel estimator, it is 
necessary to smooth the noisy estimates using some form of 
smoothing filter, such as a lowpass filter, a LS line-fitting 
filter or an adaptive linear predictor. 

The delay block in the upper part of FIG. 3 represents the 
time taken to perform the hard-decision and complex con- 
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jugation operations, while the delay block in the lower half 
of that diagram is necessary to make the feedback loop 
realisable. Usually, the latter delay is taken to be one symbol 
interval. 

In a multipath environment, the block diagram is different 
in that it would then be necessary to perform a Rake 
combining to obtain the symbol decision b^i). Apart from 
that, the mechanism of FIG. 3 for channel estimation in a 
single-path channel is identical to that used for individual 
coefficients in a multipath channel. 

Simulation results are discussed below and demonstrate 
the performance of the proposed Kalman filter detector 
compared to the conventional matched-filter detector as well 
as to Chen and Roy's RLS discussion in D. S. Chen and S. 
Roy, "An adaptive multi-user receiver for CDMA systems", 
IEEE J. Sel. Areas Comms., vol. 12, no. 5, pp. 808-816, 
June 1994. Unless otherwise stated, the system has four 
users with relative delays of 0,1, 2 and 3 chips respectively, 
random codes are used, and the Kalman filter detector 
assumes an observation noise variance of cr^-OA 

0^=0. 1 was also tried, but the results obtained were 
essentially identical to those using o 2 =04. 

FIG. 4 shows bit error rate (BER) curves for the four users 
with a processing gain of 16, whose signals are received 
with identical power levels, using both the conventional 
detector and the proposed Kalman filter detector. The noise 
variance o 2 used in the algorithm was 0.4 and the single-user 
bound is shown on the same figure for comparison. It is quite 
evident that the detector embodying the invention performs 
much better than the conventional detector, and its perfor- 
mance continues to improve with SNR even at 10 dB, 
whereas the performance of the conventional detector is 
already beginning to saturate at that SNR level. 

FIG. 5 presents BER curves obtained under identical 
conditions except that the processing gain has been 
increased to 32, to verify that the Kalman filter detector, 
although it does not perform the well-known despreading 
operation, does show improved performance with a larger 
processing gain. 

FIG. 6 gives BER curves for the four users when the 
received signal power levels are 0, 3.5, 6 and 8 dB respec- 
tively. The processing gain is 16. The solid curves are results 
obtained using the Kalman filter, while the dashed ones 
come from the conventional detector. The circles mark the 
BER's of the first user when power control is perfect, and it 
is seen from the closeness of the circles to the first solid line 
that the detector embodying the present invention is effec- 
tively near-far resistant. 

FIG. 7 compares the RLS detector discussed in 20 with 
the conventional detector and a Kalman filter detector 
embodying the present invention. It was found through 
many simulation trials that the forgetting factor 1 in the RLS 
algorithm does not affect the BER very significantly. 

The results shown in FIG. 7 were produced using 1«0.995 
The initial value of the inverse correlation matrix P was 
chosen to be the identity matrix. The relative delays between 
signals received from the four users were 0, 3, 6 and 9 chips 
respectively. The processing gain is 16. From FIG. 7 it is 
quite apparent that the RLS detector performs rather poorly, 
in fact under-performing even the conventional detector at 
low SNR's. But as already explained in the previous section, 
this is not too surprising considering that it is severely 
suboptimal when used in an asynchronous system. 

Referring to FIGS. 8 and 9, a method for multipath 
searching is provided to solve the problems associated with 
the poor performance of sliding correlators. In the method, 
as shown schematically in FIG. 8, the following steps are 
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utilised: the symbol estimates b^i), k=l, . . . , K, and 
remodulate (or respread) each one by the corresponding 
spreading code s^n); add the respread signals together, and 
subtract the result from the received signal r(n) to form e(n); 
and use e(n) as the input to a bank of K conventional 
multipath searchers, one for each user. 

This method makes use of the interference rejection 
capability of the proposed detector to provide a cleaner 
signal to the multipath searchers. Specifically, e(n) should 
contain only thermal noise when the channel multipath 
profile has not changed substantially from the last time 
channel estimation was performed. In this case, the searcher 
signals to the detector that there is no change in the channel. 
However, when the channel profile has changed 
significantly, the input to the searcher consists of the con- 
tribution from the new, uncancelled path as well as noise, but 
little MAI because it has mostly been cancelled by the 
proposed detector. The searcher will then tell the detector 
about the new path, and the detector should then act on that 
information appropriately. 

Referring to FIG. 9, the idea of generating a residual error 
signal e(n) may also be used for improved access channel 
decoding. The P access channel(s) on the reverse link of a 
cellular system are used for initial establishment of the link 
between a new mobile user and the base station, and are 
shared by all users. Correct demodulation of the access 
channels) at the base station is vital in order that new users 
may be added to the base station's list of active mobile 
stations. But again, access channels are also affected by 
MAI. Each access channel may be treated as an additional 
user inside a multi-user detector, but this only adds to 
complexity and is an inefficient use of resources since the 
access channel is not very heavily loaded, relative to the load 
carried on the traffic channels. 

The cancellation of all active users' traffic channel signals 
from the received signal leaves behind only the access 
channel signals, which can then be demodulated with a 
higher degree of reliability than using the received signal 
itself. 

What is claimed is: 

1. A method for receiving a code-division multiple-access 
(CDMA) signal combining symbol values spread by a 
plurality of spreading codes, which method comprises the 
steps of: 

(a) converting the CDMA signal to a baseband signal; 

(b) recognising symbol boundaries in every multipath for 
every user in the baseband signal to produce a delay 
estimate; 

(c) estimating an attenuation and phase rotation impressed 
on the signal by each multipath to produce channel 
estimates; 

(d) generating new symbol estimates for every user at a 
sampling rate that is greater than a symbol rate of the 
signal; 

(e) generating an estimate of the received signal using 
said symbol estimates and channel estimates; 

(f) comparing the received signal with its estimate and 
feeding estimation error back to a means for estimating 
the transmitted symbols; 

(g) sampling the symbol estimates at the estimated sym- 
bol boundaries of each user to obtain final symbol 
estimates; and 

(h) repeating steps (b) through (g) with every new sample. 

2. A method according to claim 1, wherein the channel 
estimates are used only in the symbol demodulation. 

3. A method according to claim 1, wherein the symbol 
estimates are generated using a Kalman filter. 
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4. A method according to claim 1, wherein the delay and 
channel estimation steps are performed at an update rate 
other than the sampling rate. 

5. A method according to claim 4, wherein an update rate 
lower than the sampling rate is used. 

6. A method according to claim 1, wherein the modulation 
scheme is not BPSK. 

7. A method according to claim 1, wherein a state estimate 
vector b(n) or &(n) is re -initialised at the symbol boundaries. 

8. A method for receiving a code-division multiple-access 
(CDMA) signal combining symbol values spread by a 
plurality of spreading codes, which method comprises the 
steps of: 

(a) converting the CDMA signal to a baseband signal; 

(b) recognising symbol boundaries in every multipath for 
every user in the baseband signal to produce a delay 
estimate; 

(c) estimating an attenuation and phase rotation impressed 
on the signal by each multipath to produce channel 
estimates; 

(d) generating new symbol estimates for every user at a 
sampling rate; 

(e) generating an estimate of a received signal using said 
symbol estimates and channel estimates; 

(f) comparing the received signal with its estimate and 
feeding an estimation error back to a means for esti- 
mating the transmitted symbols; 

(g) sampling the symbol estimates at the estimated sym- 
bol boundaries of each user to obtain final symbol 
estimates; 

(h) repeating steps (b) through (g) with every new sample; 
and 

multipath searching at output of a receiver, the multipath 

searching comprising the steps of: 
respreading estimated symbols transmitted by each user to 

derive respread signals; 
summing the respread signals to derived a summed signal; 
subtracting the summed signal from the received signal to 

derive a residual error signal; and 
feeding the residual error signal to a bank of multipath 

searchers. 

9. A method according to claim 8, wherein each of the 
multipath searchers in the bank of multipath searchers 
comprises one or more sliding correlators, 

10. A method according to claim 8, comprising the step of 
access channel demodulation at the output of the receiver, 
which method comprises the steps of 

respreading the estimated symbols transmitted by each 
user to derive respread signals; 
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summing the respread signals to derive a summed signal; 
subtracting the summed signal from the received signal to 

derive a residual error signal; and 
then correlating with the spreading code used on the 

access channel. 

11. A code division multiple access signal receiver com- 
prising: 

means to receive a CDMA signal combining symbol 
values spread by a plurality of spreading codes; 

means to convert the CDMA signal to a baseband signal; 

means to recognise symbol boundaries in every multipath 
for every user in the baseband signal to produce a delay 
estimate; 

means to estimate an attenuation and phase rotation 

impressed on the signal by each multipath to produce 

channel estimates; 
means to generate new symbol estimates for every user at 

a sampling rate greater than a symbol rate of the signal; 
means to generate an estimate of the received signal using 

said symbol estimates and channel estimates; 

means to compare the received signal with its estimate 
and feed an estimation error back to a means for 
estimating the transmitted symbols; and 

means to sample the symbol estimates at the estimated 
symbol boundaries of each user to obtain a final symbol 
estimate. 

12. A method according to claim 1, comprising the step of 
access channel demodulation at the output of the receiver, 
which method comprises: 

respreading estimated symbols transmitted by each user to 

derive respread signals; 
summing the respread signals the derive a summed signal; 
subtracting the summed signal from the received signal to 

derive a residual error signal; and 
correlating with the spreading code used on the access 

channel. 

13. The method of claim 1 wherein the new symbol 
estimates are successively improved. 

14. The method of claim 1 wherein the new symbol 
estimates are successively closer to the transmitted symbols. 

15. The method of claim 8 wherein the sampling rate is 
greater than a symbol rate of the signal. 

16. The method of claim 9 wherein the sampling rate is 
greater than a symbol rate of the signal. 

17. The method of claim 10 wherein the sampling rate is 
greater than a symbol rate of the signal. 
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Title page. 

Item [56], under OTHER PUBLICATIONS, "Minimization Probability of Error" should 
read - Minimum Probability of Error --. 

Column 3, 

Line 64, "code-division division" should read -- code-division --. 
Column 4. 

Line 12, delete "that is greater than a symbol rate of the signal". 
Line 35, "feed a estimation" should read - feed an estimation --. 

Column 6, 

Line 10, 'W should read ~ T ^ -. 
Line 23, "(n)] 1 " should read - (n)] T --. 
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Line 8, "K(n):=P(n)S T (n)/S(n) /^n^Ho 21 " should read 
- K(n):=P(n)S T (n)/[S(n) / 5 (n>S r (/i)+a 2 ] ~. 

Column 9, 

Line 52, "factor 1" should read ~ factor X --. 
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Line 39, "derived" should read - derive --. 
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